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Abstract 

The origin of the matter-antimatter asymmetry of the universe is one of the 
major unsolved problems in cosmology and particle physics. In this paper, we 
investigate the recently proposed possibility that split fermion models - extra di- 
mensional models where the standard model fermions are localized to different 
points around the extra dimension - could provide a means to generate this asym- 
metry during the phase transition of the localizing scalars. After setting up the 
scenario that we consider, we use a single fermion toy model to estimate the re- 
flection coefficients for scattering off the phase boundary using a more realistic 
scalar profile than previous work resulting in a different Kaluza Klein spectrum. 
The value we calculate forns/s is consistent with the mechanism being the source 
of the baryon asymmetry of our universe provided the B — L violating processes 
have an efficiency of order 10~^. 



1 Introduction 

It is clear, at least in our neighborhood of the universe, that matter is far more common 
than antimatter. Observations seem to rule out the possibility that the universe is made 
of separate regions of matter and antimatter domination due to the non-detection of 
annihilation photons from the boundaries between the regions [1]. Further, there is no 
known mechanism which would lead to such an ordering. It thus appears that there is 
far more matter than antimatter in the universe, or equivalently that the universe has 
a net baryon number, B. The origin of this baryon asymmetry of the universe (BAU) 
is one of the major unanswered questions facing cosmology and particle physics. The 
standard way to parameterize the size of the asymmetry is as the ratio of the universe's 
baryon number density to its entropy density: 

ub/s'^ 10-^°. (1) 

Given this, we require a theory which can generate a baryon asymmetry from an 
initially baryon symmetric universe. Sakharov [2] showed there were three conditions 
necessary for a theory to create such an asymmetry: 1) Baryon number violation, 2) 
C and CP violation, and 3) Departure from thermal equilibrium'. It is well known 
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'it should be noted that none of Sakharov's conditions are strictly necessary for baryogenesis to occur 
and models have been developed where each can be circumvented [3] . However, such models have thus far 
had Hmited success and at least for the type of models we consider here, all three conditions are necessaiy. 
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that these three requirements may be met by the Standard Model (SM) during the 

Electroweak Phase Transition (EWPT) via a mechanism known as electroweak (EW) 
baryogenesis [4] but unfortunately the parameters of the SM - in particular the values 
of parameters in the CKM matrix and the mass of the Higgs Particle - conspire to make 
the resulting asymmetry many orders of magnitude too small to explain that which we 
observe. However, it should be noted that the basic mechanism of EW baryogenesis - 
the CP non-invariant scattering of fermions off the boundary between two phases dur- 
ing a phase transition - remains a feasible means to generate the asymmetry if we can 
find an alternative context in which it might take place. As we will show it is possible 
that models with split fermions may provide such an alternative context. 

Split fermion models are extra dimensional models in which the SM fermions are 
confined to a brane by one or more localizing scalar fields which have a non-trivial 
vacuum expectation value (VEV) along the extra dimension. (Throughout we will 
assume, for simplicity, that there is only a single extra dimension.) However, unlike 
the standard scenario the position along the extra dimension to which the fields are 
localized varies for the different quark and lepton flavors leading to a so-called fat brane 
structure. Such models are attractive since they can provide a natural explanation for 
both the hierarchy in the SM Yukawa couplings and the suppression of proton decay 
with both these features of the SM resulting from the exponentiaUy small overlap of 
wavefunctions of fermions localized to different points. 

As was noted in [5,6], such a scenario also provides an arena where we might hope 
to meet all Sakharov's conditions and where baryogenesis may have occurred. This 
would have taken place during the phase transition where the localizing scalar fields 
acquire VEVs. In analogy to the EW baryogenesis model, during this period there are 
regions of both broken and unbroken phase. In the unbroken phase, where the scalars 
have no VEV, the SM fermions are unlocalized and the theory is fully 5d. In this 
phase the fermion wavefunction overlaps are large and we expect Yukawa couplings 
to naturally be 0(1) and baryon number violating processes to proceed rapidly. In the 
broken phase, the SM fermions are localized and the effective theory is the 4d SM - 
Yukawa couplings are hierarchical and baryon number violating processes (with the 
exception of sphaleron processes) are exponentially suppressed. Thus both B and CP 
violation can be large in the broken phase and we may have a departure from thermal 
equilibrium provided the phase transition is first order. 

To test if this model is feasible involves calculating the CP non-invariant scattering 
of SM fermions off the phase boundary and determining the size of the asymmetry in 
the reflection of particles and anti-particles. However, as is generally the case when 
we consider flelds propagating in compact extra dimensions, the problem is compU- 
cated by the presence of Kaluza-Klein (KK) modes of the SM particles. Further, the 
additional scalars we must add to the theory in order to localize the fermions introduce 
numerous, relatively unconstrained parameters. In this paper we therefore assume that 
the asymmetry generated in a realistic three generation model may be approximated 
by a calculation in a simplified toy model with just a single fermion. This approach 
has previously been considered to estimate the BAU generated in split fermion mod- 
els [5, 6] with results consistent with this mechanism being the source of the BAU. 
However, these papers consider a simpUfied profile for the localizing scalar leading to 
a markedly different KK spectrum in the broken phase than that found in a detailed 
study of this spectrum with a more reaUstic scalar profile [7]. Since the reflection 
asynmietries depend crucially on the KK spectrum in the two phases (including off- 
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shell modes), this has the potential to dramatically alter these previous results.^ 

The outline of this paper is as follows: In Section|2]we give a brief review of elec- 
troweak baryogenesis while in Section|3lwe introduce the original split fermion setup. 
We outline the details of our implementation of split fermion baryogenesis in Section 
l^and describe the simplified single fermion model and the KK modes it will contain in 
Section|5l(with some necessary results from a model with an infinite extra dimension 
contained in AppendixIXt. In Section|6lwe calculate the reflection coefficients for this 
model and find the resulting asymmetry before concluding in Section^ 

2 Electroweak Baryogenesis 

In the following paragraphs, we give a very brief outline of the electroweak baryogen- 
esis process, largely as an analogy to the split fermion baryogenesis mechanism which 
we will consider for the remainder of the paper. For a more thorough review see [4] 
and references contained therein. 

EW baryogenesis attempts to generate the BAU using only SM physics. All three of 
Sakharov's conditions may be met: 1) Baryon number is violated by non-perturbative 
sphaleron processes [10], 2) CP violation occurs due to a non-zero phase in the CKM 
matrix, and 3) there may be a departure from thermal equilibrium provided the elec- 
troweak phase transition is first order. Above the temperature of the EWPT, Tewpt, 
the Higgs particle has zero VEV and sphaleron processes may proceed rapidly. Once 
the universe cools to Tewpt, regions of broken phase, in which the Higgs has non-zero 
VEV and sphaleron processes are exponentially suppressed, begin to appear As the 
universe cools further these bubbles of broken phase expand and envelop surrounding 
regions until they fill the visible universe. As the bubbles expand, quarks from the 
unbroken phase scatter off their walls. Due to CP non-invariance, this scattering is dif- 
ferent for quarks and antiquarks and their reflection coefficients, rI^, (represent- 
ing the scattering of a left-handed quark of flavor i into a right-handed quark of flavor 
/, and its CP-conjugate process, respectively), will differ As a result, a net baryon 
number will build up inside the bubble with an equal and opposite number outside. Re- 
call, however, that B violation proceeds rapidly outside the bubble so this asymmetry 
is quickly wiped away while inside the bubble B violation is suppressed, allowing the 
asymmetry here to be preserved. Thus, once the bubbles expand to encompass the vis- 
ible universe, we may be left with the required matter-antimatter asymmetry. Ignoring 
the effect of particles from the broken phase scattering back into the unbroken phase, 
the resultant asymmetry is given by: 



where A(£') = Tr(i?^j^_RRL — /J^r^lr), and nj^^ the Fermi-Dirac distributions in 
the unbroken phase boosted to the wall frame. Taking the wall speed to be « 0.1, 
we may expand this expression to first order in obtaining [11]: 



where no{E) is the rest frame Fermi-Dirac distribution function and pr.l are the par- 
ticle momenta. 

^ [8] and [9] also use split fermion models to generate the BAU but via very different mechanisms to that 
we consider here. 




(2) 




(3) 
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Hence to calculate ns/s, we need to find values for the reflection coefficients. This 
is done by writing down wavefunctions in each phase for particles scattered off the bub- 
ble wall in terms of reflection and transmission amplitudes. If we assume the wall is 
thin (i.e. 6^ < v^Ti where 6uj is the wall thickness and r, are the timescales of pro- 
cesses taking place), continuity requires the wavefunction in the two phases to match 
on the wall. Enforcing this requirement allows us to find the reflection/transmission 
amplitudes and in turn the reflection coefficients and the resultant BAU. 

While this model of baryogenesis is attractive in that it does not require any new 
physics (beyond the as yet undetected Higgs particle), it is now ruled out, at least in 
its simplest formulation. CP violation in the SM is a small effect as parameterized by 
the Jarlskog parameter, J ^ 10"^ [12] and there has been disagreement as to whether 
it is sufficient to generate the required asymmetry once quantum effects are fully in- 
cluded [11, 13, 14]. More categorically, the Higgs mass is now restricted to be greater 
than 114 GeV [12]. Such a large Higgs mass leads to a second order EWPT [15, 16] 
preventing the required departure from thermal equilibrium. Further, with a large Higgs 
mass, the suppression of sphaleron processes below Tewpt is reduced and these pro- 
cesses can proceed at a rate sufficient to wipe out any asymmetry that is generated. 
We emphasize that these are problems related to the parameters of the SM and are 
not intrinsic to the mechanism described. It is therefore conceivable that the asymme- 
try is generated during a phase transition where some scalar other than the SM Higgs 
acquires a VEV. It is this possibility we explore in the context of split fermions. 

3 Arkani-Hamed Schmaltz Model 

In their original proposal [17], Arkani-Hamed and Schmaltz (AS) considered a fermion 
and a scalar in an uncompactified 4H-ld space with a Lagrangian: 

C = '4){x,xz){i^^'d^, + 17''' ^5 + (l){x,xz) + ■m)%l;{x,xz), (4) 

where x labels our normal 3H-ld spacetime and 0:5, the extra dimension. The profile of 
the scalar field's expectation value along the extra dimension is that of a domain wall 
centred at 2:5 — 0. Initially m is set to zero. For the case where the scalar profile is 
approximated by {(i){xc,)) — 211^x5 close to X5 = 0, the Dirac equation for the fermion 
is found to have a chiral zero mode solution: 

^1/2 

i'Lfii^s) = i^^/2y/^ exp(-^2a;^), (5) 

while the corresponding right-handed solution is not normalisable. This solution is a 
Gaussian centred at X5 — 0. If the fermion is now given a bare mass i.e. m 7^ 0, 
we find the same left-handed chiral solution now centred about x^ = — m/2/i (Figure 
Q. Thus if we generalize to the many fermion case we see we can obtain chiral ferm- 
ions localized at different points along the extra dimension if we give these fermions 
different bare masses. 

AS note two attractive properties of such a setup. Firstly, we may consider SM 
fermions where the bare masses of the SU(2) singlets and doublets differ, leading to 
their separation by a distance r along the extra dimension. Adding a Higgs to the 
Lagrangian, if the fundamental 5d Yukawa coupling between the fermions and the 
Higgs is considered to be C(l), the effective 4d Yukawa couplings are found to be 
suppressed by the small overlap between the wavefunctions for the SU(2) singlet and 
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Figure 1 : When it acquires a VEV, boundary conditions force thie scalar to take a domain wall 
or kink configuration along the extra dimension. In this environment the fermion solution is a 
'Gaussian' localized around the zero of the scalar. By giving a fermion a bare mass m, we can 
shift its wavefunction to be centred about X5 — — Giving the different fermions different 
bare masses allows us to split up the fermions along the extra dimension. 

2 2/9 

doublet of a given particle. As a function of r this suppression is found to be e^^ ' . 
Thus by considering fermions with a small range of bare masses we can naturally obtain 
the observed hierarchy in SM Yukawa couplings. 

A somewhat more involved analysis shows that a similar separation between quark 
and lepton fields can lead to the almost complete suppression of dangerous operators 
in the 5d Lagrangian (such as qqql terms) which would otherwise lead to proton decay 
occurring at unacceptable rates. 

Since the work of AS, there have been numerous variations on their original setup 
and investigations into these models' phenomenology, experimental signatures and 
consistency with the SM [5,6,8,9, 18-28]. Most importantly, models where the ex- 
tra dimension is compactified as a circle of radius R have been considered [18, 19,24]. 
Such models require orbifold boundary conditions (OBCs) in order to reproduce chiral 
fermions and an alternative means to control the overlap between the fermions since 
bare mass terms are forbidden by the OBCs. Such means include the introduction 
of non-renormalizable [18] or boundary [19] terms, non-uniform scalar-fermion cou- 
plings [19] or the addition of a second localizing scalar [24]. We will consider the last 
of these options since it also has the potential to provide the source of CP violation 
which we require (see Section|4j. 

In the context of baryogenesis, this model is particularly suggestive. It contains 
a scalar field with a non-zero VEV meaning there may have been a first order phase 
transition when this was acquired. This could potentially provide us with the necessary 
departure from thermal equilibrium. Prior to the transition, the SM fermions are un- 
localized and there are no suppressions caused by small wavefunction overlaps. Thus 
B violation may proceed rapidly and, due to the large Yukawa couplings, CP violation 
may be large. It is therefore possible that all of Sakharov's conditions are met and we 
may be able to generate the BAU providing certain cosmological conditions are met as 
we discuss in the following section. 

4 Split Fermion Baryogenesis 

In order to satisfy Sakharov's conditions and generate the BAU, our spUt fermion model 
must contain CP violation. It has recently been proposed that it may not be necessary 
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to include any new sources of CP violation and that the CKM phase could be all that 

is required to generate the BAU using split fermions [6]. We do not consider this 
possibility here and instead follow [5] and introduce a second localizing real scalar to 
our Lagrangian: 

c- ~ w^<^^ + f'i<^i)'i>3 + (6) 

In a model with three fermion generations, the will contain one physical complex 

phase, if, which will lead to CP violation. Since we anticipate all the l/^-'jl to be of 
order one, we expect the CP violation to be large. As mentioned above, the addition of 
a second localizing scalar would also allow for the fermions to be localized to different 
points in the bulk, and it has been shown [24] that it is possible to reproduce the SM 
Yukawa couplings, although the choice taken there is not the only possibiUty and the 
/s are not uniquely determined. 

Having set out the particle content of the model, the most important factor in de- 
termining its cosmological context is the temperature, Tbau> at which the locaUzing 
scalars acquire VEVs (for simplicity we assume this happens simultaneously for both 
(/)i,2) leading to a phase ttansition and the possible generation of the baryon asym- 
metry. General arguments suggest that to be consistent it is Ukely we require that 
Tbau ^ 1 /tt-R [5]. Depending on the particular split fermion model considered and the 
regions of the parameter space explored, studies searching for lower bounds on 
(resulting from the need to suppress additional flavor changing neutral currents) find 
results ranging from a few TeV to a few thousand TeV [25, 29, 30]. For our purposes, 
it is important that this suggests Tbau ^ ^ewpt and the Higgs has zero VEV during 
the epoch when we are considering baryogenesis to have taken place. This means we 
are anticipating that the BAU will be generated well above the EWPT during an epoch 
when sphaleron processes which violate B + L were proceeding rapidly. As a result, 
if only B + L were violated in our construction, any BAU generated at Tbau would 
have been completely wiped out by the time sphaleron processes were suppressed at 
the EWPT. It is therefore necessary that the asymmetry is produced by B — L violation. 
The resulting asymmetry will clearly depend on the efficiency of this process, cb-l- 

The leading order baryon violating processes allowed by the SM gauge symmetries, 
from terms of the form -^QQQL, where is the cut-off, conserve B — L and are 
therefore not relevant in this context. The lowest order processes which can violate 
B — L have the symbolic form ■^{QDmQ){QT^ L) where Dm is the covariant 

derivative, are the gamma matrices in 5D and Ad = 0, 1, 2, 3, 5. In [6] it was 
shown that assuming the theory is strongly coupled at M», naive dimensional analysis 
(NDA) leads us to estimate that these processes proceed at a rate: 

In order that the asymmetry generated inside the bubble by the CP non-invariant scat- 
tering off the wall is preserved, the opposite asymmetry outside the bubble wall must 
be washed-out hy B — L violating processes before the bubble expands to encompass 
that region of space. Using Tbau ~ l/7ri?, the efficiency of this process is therefore 
given by: 

1 / 1 \^ 

where tw ~ 4/(rti^) is the timescale associated with the thin bubble wall [6, 31] and 
there is a factor of 2 since each interaction violates B — Lhy2 units. We are unable 
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to constrain this inefficiency any further since the split fermion mechanism effectively 
blinds us at low energies to the physics at energies above Tbau and in particular the 
value taken by AI^, . It is possible that this can be better constrained by cosmological 
considerations but such an analysis is well beyond the scope of this paper 

The setup is now analogous to that of EW baryogenesis described above. Again 
the BAU will be generated by asymmetries in the scattering of quarks and antiquarks 
at the boundary between the regions where the scalars have zero and non-zero VEVs 
respectively. By analogy with Equation (|2ji, we can write: 



m.n 



(9) 

Note that since we are working in an extra dimension we must sum over not only the 
flavors of the incoming and outgoing modes, but also their KK numbers, n, m. This 
means 

A/'"'^<- = (<'''""'^<^'""' -<^""'^<^'""') - |<'^'""Vsin,, 

(10) 

and the asymmetry is difficult to calculate, even numerically. Aside from calculational 
difficulties, there is yet to be a dynamical realization of the split fermion idea which re- 
produces the Yukawa couplings of the SM while suppressing proton decay. As a result, 
even if we could perform the calculation, we do not have a fully realistic model to pe- 
form it on. We therefore simplify the problem by calculating the reflection coefficients 
in a model containing only a single fermion flavor. In doing this, we anticipate that 
they will carry over to the three generation case with corrections of 0(1) [6]. Having 
made this simplification, the equation for the asymmetry reduces to: 

1 ^ 10-3 ^ ^ dE^^^E^^^ _ A«™(£;), (H) 



Tbau J Stt Tbau 



with A"™(£') ^ |7^"™p (where 7^"™ are the single fermion reflection coefficients) 
since we anticipate an order one phase. Our problem is now reduced to finding the 
reflection coefficients. To do this we must find the fermion wavefunctions for all KK 
modes, in both phases. In this single fermion framework, neither of the qualitative 
effects of adding a second localizing scalar - a non-removable complex phase and 
fermion splitting - is present. We are therefore able to further simplify the calculation 
by finding the coefficients in a model with only one scalar We carry out this calculation 
in the follow sections. 



5 A Single Fermion Toy model 

In our toy model, we have a single massless fermion and a single real scalar moving in 
a 4h-1 dimensional spacetime. The Lagrangian is 

£ = ^/j{x,X5){i-f''d^, - j^d5 - f<j){x,x^))^lj{x,x^) + ^5''0(x, X5)5,,(/)(x, 2:5) 
- i9^0(a;,a;5)550(x,X5) - ^(^^(x, xg) - i;')'. (12) 
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In order to obtain the chiral fermions, which we require to reproduce the SM, we 
compactify the extra dimension as an orbifold /7i2, which imposes the boundary 
conditions: 

(j){x, — ^5) = — 0(a;, xs), </)(a;, ttR + x^) = — —-kR + x^), 
ipix^-xs) ^ -f5i>{x,X5), iI:{x,ttR + x^) ^ -i^ili[x,-TiR + x^), (13) 

where R is the radius of the extra dimension. 

As with any theory where particles can propagate in compact extra dimensions, the 
effective 4d Lagrangian will contain a tower of KK modes: 

'4>{X,X^) = ^ {i>nR{x)inn{Xh) + tpnAx)S.nL{Xb)) , (14) 
n=0 

where we have also performed a chiral decomposition. This leads to a Dirac equation 
for the extra dimensional part of the wavefunction: 

{-d5 + f{(j){xr,)))£_,n^ = m^jnR- (15) 

Above the breaking scale, where {(j){x5)) — 0, these equations are easily solved 
and with the OBCs il3\ just give the usual KK chiral modes which we label Xn- 

''"^hn-^x„f'''""^(S V rn,.-l neN. (16) 
Xnx,(a;5) J ^" \ smnx^/R ^" R' 

At a temperature Tbau ~ l/nR the scalar may undergo a phase transition and 
acquire a VEV. It is at the boundary wall between these two phases that we envisage 
the BAU being generated. The OBCs, Eq. (I13> . requiring (j) to be odd about the fixed 
points, lead to a double kink solution along the extra dimension, given approximately 
by [18]: 



(0(a;5))=«tanh|^y^a;5jtanh|^y^(^i?-X5)j + ©(e-'^^^). (17) 

Clearly this approximation is valid so long as the brane is much narrower the size of 
the extra dimension^ i.e. 

-r^ < ttR. (18) 

We note that this is the same condition required for a split fermion model to effectively 
suppress Higgs Yukawa couplings and proton decay. In [5, 6] this scalar profile was 
approximated by a step function. 

Combining equations il5\ and substituting in Eq. illl for the scalar profile we 
obtain two second order differential equations [7]: 

{-dl + VL)^rn,^m^U,, (19) 



^An exact solution has been found for the scalar VEV profile (as well as its KK modes and those of the 
fermion) [28] and as expected it is in very good agreement with the approximation in this limit. 
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where, 



f tanhuinR — X5) tanhwxs \ 

Vb,l = Tuw , 2 / n ^ 

V cosh UX5 cosh u{TrK — X5)/ 

+ w'^ tanh^ ux^tanh^ u{ttR — X5), (20) 



and. 




w = fv, u = d —. (21) 

In general, we can't solve these equations but in the narrow brane case the solu- 
tions can be related to the solutions about a scalar kink in an infinite extra dimension. 
These solutions were derived in [7, 18] and we reproduce them in AppendixlAl In the 
following paragraphs we regularly refer to objects defined therein. 
Note, in the narrow brane limit: 

, , I Vl^^\x5) near 2:5 = 0, 

Vr{x^) ~ { /J^)^ 

ttR) near 2:5 — nn, 

\r ( 'L v-s; nearxs^O, 

ttR) near 2:5 — nR, 

thus we can approximate the solutions by those for the infinite case, Eqs. (|37},(|40l. 
As in the infinite case we find both bound solutions localized around the fixed points 
and unbound solutions which propagate throughout the extra dimension. However, the 
addition of OBCs further restricts the allowed solutions by requiring the right-handed 
(left-handed) solutions to be even (odd) about both fixed points. 

Thus, for the bound modes, when n is even, ^^^(xs) ~ CnnH^s) around 2:5 = 
while the solution around 2:5 — ttR is suppressed since it would be odd i.e. (x^) — 
about 2:5 — ttR. Similarly, when n is odd, the right-handed solution is suppressed 
around 2:5 = but takes the form Cnnixs) ~ £,ii^\x5 — nR) about 2:5 = nR. The 
matching of the allowed solution to zero at the mid-point between the branes will be 
unproblematic since the solution will have already asymptoted to zero in the narrow 
brane context we are considering. Similar arguments apply for the left-handed solution 
with there again being solutions localized to 2:5 = (2:5 = irR) for n even (odd) 
(although when n — there is no allowed solution, ensuring we get the chiral zero 
mode we require). This is summarized graphically for the first few modes in Figure 
|2] The expression for the masses of these modes carries straight over from the infinite 
case: 

m? = 2nuw ~ n^v^ , whereO<n<— . (23) 

^" u 



Note that provided w > u (i.e. / > y ■^), there will be massive states localized at the 
fixed points. Such modes were not present in the analysis of [5, 6] as their simplifica- 
tion of considering the scalar kink as a step function is equivalent to taking A ^ 00 
thereby removing these localized massive modes. Since the reflection asymmetries are 
essentially determined by the KK spectra in the two phases, this difference is poten- 
tially significant (especially since these modes are the lightest KK modes in the broken 
phase). 

For the unbound modes, there exist both odd and even solutions of each chirality 
in the infinite extra dimension case (see Eq. gO}) so we can always find a solution 
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n=0 




Figure 2: The bound modes, m < w, are localized about the fixed points of the extra dimension. 
The even KK modes are found about x-, = while the odd modes are found at xr, = nR. OBCs 
require that the right-handed modes (solid lines) be even functions about the fixed points, while 
the left-handed modes (dotted lines) must be odd. There is no normalisable left-handed zero 
mode ensuring the SM derived from this model will be chiral as required. 



satisfying the OBCs about each fixed point: 

r \ - ) CSeven(a;5), < X5 < TtR/2, 

^mR^^S) — S (00) . I^^ D/o ^ / D 

[ Cmt,even(.^;5 - TT-K), 'kR/2 < < TtR, 

{X,) - I _ ^^^^ ^^^^ ^^^^ (25) 

Continuity of the wavefunction means we must match these solutions at the mid-point 
between the fixed points, = itR/2. For the right- and left-handed solutions respec- 
tively, this gives: 

d^!eve„M/2) = C^^leve„(-^^/2), (26) 

CSodd(-^/2)-CSodd(--^'/2)- (27) 

In order that such a matching can be done in general, the phases of the two wavefunc- 
tions must match, which can always be achieved by an appropriate choice of normal- 
ization constants. The matching only occurs for certain values of m,^^^ , leading to a 
discretized tower of allowed unbound modes. Note these matching conditions differ 
from those found in [7]. In general we cannot solve these equations analytically and 
we must resort to numerical techniques. Having done so, we will have an infinite tower 
of modes of each chiraUty which we call Cur^l ■ 



6 Reflection Coefficients and Baryon Asymmetry 

Having found the solutions to the Dirac equations in both the broken and unbroken 
phases, we now proceed to calculate the reflection amplitudes for the scattering of 
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incoming KK fermions off the phase boundary. This involves matching the wavefunc- 
tions of the fermions in the broken and unbroken phases on the boundary between 
them, in exact analogy with the procedure for EW baryogenesis described in Section 
121 As in that case, we may Lorentz transform away any components of the momentum 
of the incoming particle along uncompactified dimensions parallel to the wall (that is, 
along the xi and X2 directions if we define the xs-axis to be normal to the wall) and 
position the wall at X3 — 0. We cannot, however, transform away momenta along the 
extra dimension since here Lorentz invariance is broken by the orbifold fixed points. 
As a result we are matching the wavefunctions not only at a point on the boundary but 
right around the extra dimension. 

We consider in turn an incoming particle of each KK mode, Xn', with helicity 
h = ±1. This can be reflected back in to the unbroken phase as Xn with amplitude r„, 
or transmitted to the broken phase as or ^„ with amplitudes t(^^ and t^^ respectively. 
The wavefunction matching condition on the bubble wall is then: 



, n=0 



00 

E 

n=l 



hk^-'x 



hkl^'x 




iEt 



X3=0 



E 



hktCnn{x5) 



E 



E - hk^zini^ix'i) 



-iEt 



(28) 



where \y\ is the largest integer smaller than y, E is the energy of the incoming particle 

2 are the particle momenta in the 2:3 -direction 



and fcg"'^"''^ 



where we recall that the masses must be found numerically for the (^„s. 

To solve this equation, we find Fourier series expansions for the ^„ and Cn. We may 
then equate the Fourier coefficients on the two sides term by term. Since this must be 
done numerically, it is necessary to truncate the infinite sums in Eq. ( I28> at and N'^ 
respectively. For sufficiently high truncation, this leaves us with an over-determined set 
of simultaneous equations in the rs and ts which we may solve by the pseudo-inverse 
method. 

We can then find the reflection and transmission coefficients: 



(29) 




Unitarity should ensure that 

7e"' [E) 



ATX 



N 



[E) = ^Ti" 

n=0 



\w j u\ 
n—1 n— 



= 1, 



(30) 



(31) 



(32) 
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for all n' and E and this is indeed the case provided we choose N^-'^ sufficiently large. 
Plots of 7?."" as a function of the incoming particle's energy are shown in Figure |3l 
for two choices of values for the parameters /, A, v. In these plots the effect on the 
reflection coefficients of additional modes coming on-shell as the incoming particle's 
energy is increased is clearly visible. It is also apparent that the choice of parameters 
has a significant impact on the values the various reflection coefficients take, suggesting 
the value we find for ub/ s may be strongly parameter dependent. 




I ^ ^ — ' ' ' ' ' ' 1 

2 4 6 8 10 12 14 16 

ER 

1 I 1 1 1 1 1 1 1 1 

0.9 - /""n 

0.8 - / i 

0.7 - ; i 



0.6 - 




Figure 3: Plots of T?,^'" (solid line), T?.^'" (dashed line) and Ti?'^ (dotted line) as a function of the 
incoming particle's energy for n = 10/7?,/ = 1, A = 1 (top plot) and w = 20/R,f = 1,A = .1 
(bottom plot). 

However, when we use Eq. il 1\ to calculate the BAU we find ub/s ^ lO^^eB-L, 
relatively independently of the choice of values for the parameters, provided we remain 
in the narrow brane regime of Eq. dlSl l. Further, as long as Tbau is of the same order 
as 1/ ttR, its precise value does not change the order of magnitude of the resulting 
asymmetry. This number is a couple orders of magnitude smaller than that found in 
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similar calculations [5, 6] where the scalar kink is taken as a step function and there are 
no localized massive modes in the broken phase'*, confirming that these extra modes 
have a significant effect on the result. However, it is of the same order of magnitude as 
another calculation of rt^ /s in [6] . This calculation considers a three generation model, 
although it does not implement a quark-lepton separation to suppress proton decay. It 
uses a step function as the scalar profile and uses a different method to calculate the 
reflection asymmetries, exploring a different region of the parameter space, namely 
where v <^ 1/R. Interestingly, a small v would require us to take A large to remain in 
the narrow brane regime. Recall that taking A to be large is precisely the limit in which 
our choice of scalar profile tends towards a step function meaning the choice of scalar 
profile may not be as crucial in this region of parameter space. 

In order that this mechanism generate the known value of ub/s ^ 10^^" requires 
es-L ~ 10^^ which by Eq. (|8} gives 

M, ~ ^. (33) 

Given the relative simplicity of our construction and the NDA used to derive the ex- 
pression for es-L, one should perhaps not take this relation too seriously. However, it 
does suggest that in order for split fermion baryogenesis to be viable the fundamental 
energy scale. A/*, at which B — L violating theory applies can not be orders of magni- 
tude higher than the scale, Tbau, at which fermion localization occurs and the BAU is 
generated. 

The major uncertainty in our calculation is the degree to which our single fermion 
result carries over to a realistic three generation model. The issues here can be sepa- 
rated into two types. Firstly, there are those fundamental to going to a three generation 
model - to what extent does mixing between different fermion flavors change the re- 
flection coefficients? On this front it is suggestive that [6] find a value for ub/s in 
a three generation calculation which agrees with our result, albeit in a very different 
region of parameter space. The second issue is not fundamental to a three generation 
model but results from what else we are requiring the spliting of the fermions to pro- 
vide - a suppression of proton decay and an origin for the hierarchy in Higgs Yukawa 
couplings. As discussed early, this requires the various fermions to be localized to 
different points in the extra dimensional bulk, with quarks separated from leptons (to 
suppress proton decay) and SU (2) l doublets from singlets (to provide a Yukawa hier- 
archy). Having fermions not only localized at the fixed points but also in the bulk may 
have some important effects. Unfortunately, there has yet to be a realistic split fermion 
model which dynamical localizes the fermions in such a way as to achieve both these 
outcomes. It may be that some additional physics is required to achieve this and until 
we have such a dynamical realization of the idea, we are unable to cosider how well 
our toy model result generalizes to it. 



7 Conclusion 

Split fermion models are attractive extra dimensional constructions which have the 
potential to explain both the hierarchy in SM fermion masses and the suppression of 
proton decay. In this paper, we have carried out a calculation to estimate the baryon 

''Although [6] assumes the necessary CP violation comes from the CKM phase of the SM, while we 
assume it is from the extended scalar sector, a compaiison of results is still possible since in both cases the 
violation is 0(1). 



13 



asymmetry generated in such models during the phase transition where the locaHzing 
scalar acquires a VEV. We have improved on previous similar calculations by consid- 
ering a more realistic localizing scalar profile and resultant KK spectrum. For a range 
of parameters and transition temperatures, we find n^/s ^ 10~^eB_L compared with 
a known value of ub/ s ~ 10"^*'. This value is a couple orders of magnitude smaller 
than the past calculations but still consistent with this mechanism being the source of 
the BAU provided the scale of B — L violating physics is not too much higher than 
that of the fermion localizing physics. We have made numerous simplifications the 
most significant of which is calculating the reflection coefficients in a single fermion 
framework and anticipating that the reflection coefficients found will carry over to three 
generation models. It is interesting that our result for the generated BAU is of the same 
order as a three generation perturbative calculation [6] but this calculation is done in a 
very different region of parameter space and does not separate quarks and leptons as 
a realistic model must. While our results are suggestive, clearly more work is neces- 
sary - both to construct more realistic dynamical implementations of the original split 
fermion model and in refining our calculations to better include all the relevant physics 
- in order to draw any firm conclusions as to whether this mechanism provides a viable 
explanation of the matter-antimatter asymmetry of the universe. 
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A Split Fermion Kaluza-Klein modes in an infinite ex- 
tra dimension 

The Kaluza-Klein modes for a fermion in the presence of a scalar kink in an infinite 
extra dimension have previously been calculated by [7, 18]. Here we quote the key 
results. 

In the case of an infinite extra dimension the scalar VEV is. 



(0(2:5)) - V tanh ^^5j , (34) 

where we now have only a single fixed point at X5 = 0. Proceeding as we did in 
Section |5] we derive second order differential equations, 

i-d! + vt^)U.=m%^., (35) 

where, 

Vj^ ]^ — ^uw 2 1" tanh ux^. (36) 

cosh UX5 

With some changes of variables, these equations can be written in the form of the 
standard hypergeometric differential equation [32] which can then be solved. There 
exist both bound (m^ < w^) and unbound (m^ > w^) solutions. The bound solutions 
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are found to be: 



4f (a;5)=Ar„ 



1 



cosh UX5 
1 

cosh uxs 



2Fi[€ ,e+- + l,e+l; 

u u 

2Fi{e + l,e+-,e+l; 



1 



X -|- g2Ma;5 
1 



(37) 



where Nma^L normahzation constants, 2-^1(0, b, c; z) is the hypergeometric func- 
tion [32] and we have defined 



(38) 



In order that these solutions be normahzable, we have the additional condition that 



= 2nuw — r?v? , where < n < — , 



(39) 



leading to a discretized tower of allowed modes. For even n, is even while is 
odd. For odd n, the reverse is true. Further, when n = 0, only the right-handed mode 
is normalizable meaning we obtain the chiral massless fermion we require if we hope 
to reproduce the SM. 

For the unbound case, there are two linearly independent solutions for each chiral- 
ity, which we may conveniently write as even and odd wavefunctions, 
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cosh UX5 



2F1 e 



w w 

-,e+- 
u u 
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cosh uxz 
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\ cosh ux^ 
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Cmi,!odd(^5) = -^mi.odd 



(2e"-^)^2i^i ( - + 1, - , -e + 1; \ 
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cosh UX5 
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where. 



r(e)r(e + i) 



-e-l- 1; 



1 



]_ _|_ g2MX5 



(40) 



r(e)r(i-e) 



T{e-w/u)T{e + w/u+l)' T{w/u + l)T{-w/u)' 

D,- r(e)r(e + i) 



— . (41) 

T{e-w/u + l)T{e + w/u) ' 

In this case, there is no restriction on m and there is a continuum of modes. As we 
would expect, as x^ — > 00, these tend towards the usual plane wave solutions. 
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